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Abstract 

The inelastic production of J / ip mesons in ep collisions has been studied with the 
ZEUS detector at HERA using an integrated luminosity of 109 pb -1 . The J/ip 
mesons were identified using the decay channel J/ip —* fi + fi~. The measurements 
were performed in the kinematic range 2 < Q 2 < 80 GeV 2 , 50 < W < 250 GeV, 
0.2 < z < 0.9 and — 1.6 < Yj a b < 1.3, where Q 2 is the virtuality of the exchanged 
photon, W is the photon-proton centre-of-mass energy, z is the fraction of the 
photon energy carried by the J / ip meson in the proton rest frame and is the 
rapidity of the J ftp in the laboratory frame. The measured cross sections are 
compared to theoretical predictions within the non-relativistic QCD framework 
including colour-singlet and colour-octet contributions, as well as to predictions 
based on the /c^-factorisation approach. Calculations of the colour-singlet pro- 
cess generally agree with the data, whereas inclusion of colour-octet terms spoils 
this agreement. 
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1 Introduction 



Inelastic production of charmonium can be described in two steps. The first step is the 
creation of a cc quark pair, a process which can be calculated in perturbative Quantum 
Chromodynamics (QCD). The second step is the formation of the J /if) bound state, which 
occurs at long distances and is described by phenomenological models. 

When charmonium production was first investigated at CERN [1] and Fermilab [2] it was 
hoped that the production cross section could be used to determine the gluon density in 
the proton, because the gluon density directly enters the cross-section calculation. This 
idea was encouraged by the qualitative agreement of the data with the predictions of LO 
QCD within the framework of the colour-singlet model (CSM) [3] in which the cc pair is 
formed in a CS state identified with the J/if>. Later data irompp collisions at Fermilab [4] 
indicated that the CSM is not able to describe J /if) production at large J /if) transverse 
momenta, and hence that there may be significant contributions from higher orders in 
QCD or from the production of cc pairs in colour-octet (CO) states, which evolve into 
J / if) mesons via radiation of soft gluons. 

Models have been developed in the framework of non-relativistic QCD (NRQCD) [5] in 
which CS and CO contributions coexist. The transition of the coloured cc pair, with 
a given angular momentum, into a J/ if) is described in terms of long-distance matrix 
elements tuned to experimental data (hadroproduction of J /if) mesons or 5-meson de- 
cays to J /if))- As well as describing high-p-r charmonium production in pp collisions [4], 
NRQCD calculations including CS and CO contributions are also consistent with the data 
on production of J/ if) mesons in 77 interactions at LEP2 [6]. However, J /if) polarisation 
data from CDF [7] are inconsistent with NRQCD predictions. Comparisons with the 
decay angular distributions measured in e + e~ collisions at BaBar [8] and Belle [9] are 
inconclusive. 

The production of J /if) mesons in ep collisions at HERA is sensitive to both CS and 
CO contributions. The CS mechanism is expected to be the dominant contribution at 
intermediate values of the inelasticity variable, z < 0.7. In the proton rest frame, z is 
the fraction of the virtual photon energy transferred to the J /if). The CO mechanism is 
expected to be dominant at high z. However, large contributions from the elastic and 
diffractive proton-dissociative J /if) production processes [10,11] are also present at z ~ 1. 

Inelastic J /if) production at HERA was studied previously in the photoproduction regime 
(photon virtuality Q 2 ~ 0) by the HI [12] and ZEUS [13] collaborations. The leading- 
order (LO) NRQCD calculations and the next-to-leading-order (NLO) CSM predictions 
are both consistent with the data. Inelastic J / if) production in the deep inelastic scattering 
(DIS) regime (Q 2 > 1 GeV 2 ) has been studied by the HI collaboration [14]. In this case, 
the LO NRQCD predictions overestimate the data, whereas the LO CSM expectations 
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underestimate them. The shape of the differential cross sections are reasonably well 
reproduced by both models, except for the inelasticity distribution in the case of LO 
NRQCD and for the distribution of the J /if) transverse-momentum squared in the photon- 
proton centre-of-mass system in the case of the LO CSM. NLO CSM predictions are not 
available in the DIS regime. 

Models in the framework of the semi-hard or /^-factorisation approach [15] are also 
available. In these models, based on non-collinear parton dynamics governed by the 
BFKL [16] or CCFM [17] evolution equations, effects of non-zero gluon transverse mo- 
mentum are taken into account. Cross sections are then calculated as the convolution of 
unintegrated (transverse-momentum dependent) gluon densities and LO off-shell matrix 
elements. These models [18-20] succeed in describing the spectra of different quarko- 
nium states at Fermilab and J/ if) meson production at HERA, as well as the quarkonium 
polarisation properties measured both at Fermilab and HERA. 

This paper presents a measurement of inelastic J /if; production in DIS and comparison 
of the data with NRQCD and models in the fc^-factorisation approach. The reaction 
ep — > eJ/ipX is studied for Q 2 > 2GeV 2 . The photon-gluon fusion process leading to 
a J /if; in the final state is assumed to be the dominant mechanism. Other contributions 
arise from the production of if)' mesons which subsequently decay to a J/if), and from 
J /if; and if)' production from the resolved photon process, where the exchanged photon 
acts as a source of partons; the latter process, however, is suppressed at high Q 2 . In 
addition, beauty production at high Q 2 with subsequent decay of a B meson into a J/if> 
also contributes to the measured cross section; this process is suppressed due to the small 
beauty cross section and the small B to J/ if; branching ratio. 

Inelastic J /if; production at large Q 2 has a smaller cross section than in photoproduction 
but presents several interesting aspects. The contribution from the CO model is expected 
to be more significant; both the CO and the CS predictions should be more accurate due 
to the higher scale in the interaction. Also, backgrounds from diffractive processes are 
reduced at high Q 2 . The measurements presented here are in a larger kinematic range 
than those previously published [14]. A measurement of the hadronic final state, X, is 
presented for the first time. 

2 Experimental set-up 

The data sample used in the analysis was collected with the ZEUS detector in the years 
1996-2000 and corresponds to an integrated luminosity £ = 108.8 ± 2.2 pb _1 . During the 
1996-97 data taking, HERA operated with positrons of energy E e = 27.5 GeV and protons 
of energy E p = 820 GeV, corresponding to a centre-of-mass energy y/s = 300 GeV (£300 = 
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32.7 ± 0.6 pb" 1 ). In the years 1998-2000, HERA collided electrons or positrons with 
protons of energy E p = 920 GeV, corresponding to \fs = 318 GeV (£318 = 76.1 ±1.6 pb" 1 ). 
The cross sections presented here were corrected to y/s = 318 GeV using the Monte Carlo 
(MC) simulation described in Section HI 

A detailed description of the ZEUS detector can be found elsewhere [21,22]. Charged 
particles are tracked in the central tracking detector (CTD) [23], which operates in a 
magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists of 
72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle 1 
region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks is 
ct(j)t)/pt = 0.0058px , ©0.0065©0.0014/p T , withp T in GeV. Energy deposits are measured 
in the high-resolution uranium-scintillator calorimeter (CAL) [24] which consists of three 
parts: the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each 
part is subdivided transversely into towers and longitudinally into one electromagnetic 
section (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAC). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a(E)/E = 0.18/a/E for electrons 
and a(E)/E = 0.35/ vE for hadrons with E in GeV. The scattered electron 2 identifica- 
tion is performed by combining information from the CAL, the small-angle rear tracking 
detector (SRTD) [25,26] and the hadron-electron separator (HES) [27]. Muon identifica- 
tion is performed by finding tracks in the barrel and rear muon chambers (BMUON and 
RMUON) [28] or minimum-ionising energy deposits in the CAL, matched to CTD tracks. 
The muon chambers are placed inside and outside a magnetised iron yoke surrounding 
the CAL. The barrel and rear inner muon chambers, used in this analysis, cover polar 
angles from 34° to 135° and from 135° to 171°, respectively. 

The luminosity was determined from the rate of the electron-proton bremsstrahlung pro- 
cess, ep — > ejp, where the photon was measured by a lead-scintillator calorimeter [29] 
located at Z = —107 m. 



1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 

2 Here and in the following, the term "electron" denotes generically both the electron (e~) and the 
positron (e + ). 
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3 Event selection and reconstruction 



3.1 Event selection 

A three-level trigger system was used to select events online [22,30]. The first- and 
second-level trigger selections were based on the identification of a scattered electron in 
the CAL, as described in detail elsewhere [26]. The third-level trigger selection required 
both a scattered electron in the CAL and a track segment reconstructed in the barrel or 
rear inner muon chambers. 

Additional requirements were imposed in the offline selection in order to suppress the 
photoproduction background and select inelastic events with a J/ip candidate identified 
by the presence of a pair of oppositely charged muons. In addition, the Z position of the 
reconstructed vertex was required to lie within 50 cm of the nominal interaction point. 

3.2 Reconstruction of DIS kinematic variables 

A scattered electron candidate, identified from the pattern of the energy deposits in the 
CAL [31], was required. The electron position measurement of the CAL was improved 
using information from the SRTD and the HES. To ensure full containment of the elec- 
tromagnetic shower, the electron impact position on the inner face of the rear calorimeter 
was required to lie outside the box \X\ < 13 cm, \Y\ < 7 cm. The energy of the scattered 
electron was required to be greater than 10 GeV. 

The photon virtuality, Q 2 , was reconstructed from the polar angle and energy of the 
scattered electron and was required to be in the range 2 < Q 2 < 80 GeV 2 . The Bjorken 
variable, y = (P ■ q)/(P ■ k), where P, q and k are the four-momenta of the incoming 
proton, exchanged photon and incoming electron, respectively, was reconstructed with 
the E method [32]. Monte Carlo studies showed this method to be the most precise in the 
selected phase space region. The photon-proton centre-of-mass energy, W, calculated 
from W 2 = ys — Q 2 , was restricted to the range 50 < W < 250 GeV. 

Conservation of energy, E, and longitudinal momentum, pz, require 5 = Ej (Ei — pz,i) — 
2E e = 55 GeV, where the sum runs over all the particles in the final state. The exper- 
imentally reconstructed quantity <5 mcas = E» (Ei — pz,i) was calculated, where the sum 
runs over all final-state energy-flow objects [33] (EFOs) which combine the information 
from calorimetry and tracking. Only events with 40 < 5 mcas < 65 GeV were kept. This 
cut reduces background from photoproduction events, where the scattered electron is not 
detected, and removes DIS events with large initial-state radiation, where the incoming 
electron radiates a high-energy photon before the interaction and the photon escapes 
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detection in the rear beam hole. To reduce background from photoproduction events fur- 
ther, the condition y e < 0.95 was applied, where y e indicates the value of y reconstructed 
from the scattered electron energy and polar angle. In order to ensure an accurate re- 
construction of the final state, it was also required that the value of y, obtained with the 
Jacquet-Blondel method [34], be larger than 0.02. 

3.3 J/ip reconstruction 

The oppositely charged muons from the J/ip decays were reconstructed in the CTD. Each 
track considered in the analysis was required to be fitted to the event vertex, to reach at 
least the third superlayer of the CTD and to have transverse momentum > 100 MeV; 
this guarantees good reconstruction quality. At least one of the tracks from the J / ip decay 
had to match a segment in the inner muon chambers and the other had to match a CAL 
cluster with an energy deposit consistent with the passage of a minimum ionising particle 
(mip). To ensure high muon-identification efficiency and purity, the track matched with 
the segment in the barrel (rear) inner muon chambers was required to have px > 1.4 GeV 
(p > 1.8 GeV, where p is the track momentum). For the track matched to the mip cluster 
in the CAL, the cut p > 1 GeV was applied. The muon identification and reconstruction 
efficiencies were estimated separately for muons reconstructed in the BMUON, RMUON 
and CAL using independent samples of dimuon events. The efficiency for tracks matched 
to the inner muon chambers varies from 35% for p? ~ 1.4 GeV to 60% at high transverse 
momentum for the barrel inner muon chambers and from 50% for p ~ 1.8 GeV to 65% at 
high momentum for the rear inner muon chambers. The efficiency for tracks matched to 
a mip in the CAL was 92%. 

The J/ip rapidity in the laboratory frame, defined as Yi a b = ^^[(E^+pz,^) / (E^ — pz,ip)], 
where and pz^ are the energy and longitudinal momentum of the J/ip meson, was 
limited to the region —1.6 < Yj a b < 1.3, where the acceptance is high. 

The inelasticity of the J/ip meson, z — (P ■ p^)/(P ■ q), where p^ is the four-momentum 
of the J/ip, was reconstructed using the expression 

_ E^p — pz^, 
2E e yx ' 

where y-% = (Ei — pz,i)/5 mcas and the sum in the numerator runs over all EFOs not 
associated with the scattered electron. According to MC studies, the average resolution 
in z is 10%. The inelasticity was restricted to the range 0.2 < z < 0.9. The lower z cut 
removes the region of high non-resonant background due to fake muons and the upper z 
cut removes elastic J / ip events and suppresses diffractive J / ip production with dissociation 
of the proton. In order to suppress further the latter background, the following cuts were 
applied: 
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• the analysis was restricted to events with an energy deposit greater than 1 GeV in a 
cone of 35° along the outgoing proton beam direction (excluding calorimeter deposits 
due to the decay muons); 

• the event was required to have at least one track in addition to those associated with 
the two muons and with the scattered electron. 

Figure H] shows the invariant mass, M^^- , distribution of all selected muon pairs. The dis- 
tribution was fitted in the intervals 2.5 < M /i + At - < 3.6 GeV and 3.8 < M M + M - < 4.5 GeV 
with a function taken to be the sum of a "modified" Gaussian, to describe the signal, and 
a linear function, to describe the non-resonant background. The range 3.6 < M M + M - < 
3.8 GeV was excluded to avoid any overestimation of the background due to the ip' state. 
The modified Gaussian function had the form: 

Gauss mod oc exp[-0.5 • x ^V(^-^\ 

where x = |(M M + M - — M )/a\. This function was introduced to take into account the non- 
Gaussian tails of the resonant signal. This functional form describes both data and MC 
signals well. The position of the Gaussian, M , the signal width, a, as well as the number 
of signal events were free parameters of the fit. The fit yielded a peak position of Mq = 
3098 ± 3 MeV, in agreement with the PDG value [35], and a width of a = 35 ± 3 MeV, 
in agreement with the MC estimation of the detector resolution. The number of J/ip 
mesons was 338 ± 25. 

4 Monte Carlo models 

Inelastic J / if) events were generated using the Epjpsi [36] MC generator. Epjpsi incorpo- 
rates the photon-gluon fusion process at LO, with initial- and final-state parton showers 
performed according to the colour-dipole model as implemented in Ariadne [37]. J/ip 
mesons were produced in the framework of the CSM. The GRV98 [38] parton distribution 
functions were used. The scales for the evaluation of the strong coupling constant and the 
proton structure function were set to the centre-of-mass energy in the 7*5 frame. The 
hadronisation was performed with the Lund string model [39]. The Epjpsi MC predic- 
tions were reweighted to the data in Q 2 and p^, where p5? is J/i/j transverse momentum 
squared in the 7]? centre-of-mass frame. The helicity parameter in the Epjpsi MC was 
set to zero; this hypothesis is supported by the data [13]. 

Signal events were also generated using the Cascade MC program [40]. Cascade in- 
corporates the off-shell matrix elements for the photon-gluon fusion process at LO. The 
initial-state parton shower is generated according to the CCFM evolution equations [17]. 
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The J ft]) mesons were produced in the framework of the CSM. The gluon density, un- 
integrated in transverse momentum, kx, was obtained from an analysis of the proton 
structure functions based on the CCFM equations [41]; in the event generation the gluon 
density used corresponds to the set named " J2003 set 2" . In Cascade the hadronisation 
was also performed with the Lund string model. 

Events with diffractive dissociation of the proton, ep — > eJ/tpN, where N is a low mass 
state with the quantum numbers of the proton, were simulated using the Epsoft MC 
generator [42], which has been tuned to describe such processes at HERA [43]. Proton- 
dissociative events were also simulated with the Diffvm [44] MC generator. Diffvm has 
a more detailed simulation of the final state than Epsoft. 

J/ip mesons originating from B-meson decay were simulated using the Rapgap MC 
generator [45], via the photon-gluon fusion process, 7* g — > bb; the beauty-quark mass 
was set to 4.75 GeV. The CTEQ5L [46] parton distribution functions were used. The B 
to J/ip branching ratio in Rapgap was set to the PDG value [35]. The MC prediction 
was normalised to the measured beauty cross section in DIS [47]. 

All generated events were passed through a full simulation of the ZEUS detector based 
on Geant 3.13 [48]. They were then subjected to the same trigger requirements and 
processed by the same reconstruction programmes as for the data. 

5 Cross-section calculation 

Prior to the cross-section calculation, the residual diffractive proton-dissociative back- 
ground was subtracted. Although such events are produced at z ~ 1 and the inelasticity 
was restricted to 0.2 < z < 0.9, some diffractive events migrate into the data sample due 
to the finite z resolution. The reconstructed track multiplicity distribution was fitted to 
the sum of inelastic (Epjpsi) and diffractive (Epsoft) MC predictions. The fit yielded 
a contribution of 6 ± 1% from proton dissociation for the whole sample. The proton- 
dissociative contributions were subtracted bin-by-bin from all measured cross sections 
according to the Epsoft predictions normalised to the above fraction. 

The number of J/ip mesons reconstructed in the kinematic region 2 < Q 2 < 80 GeV 2 , 
50 < W < 250 GeV, 0.2 < z < 0.9 and -1.6 < Y Uh < 1.3, after subtraction of the proton- 
dissociative admixture, was compared to the predictions of the Epjpsi MC generator. 
The results are shown in Fig. |2] for z, Q 2 , W, p^, the J/ip rapidity in the jp frame 3 , Y*, 
and M x , where Mx is the invariant mass of the final state excluding the J/ip and the 
scattered electron. 

3 In the 7p centre-of-mass frame, the photon direction was chosen to be the "forward direction" . 
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Data were corrected bin-by-bin for geometric acceptance, detector, trigger and recon- 
struction inefficiencies, as well as for detector resolution, using the Epjpsi MC gener- 
ator. The acceptance, Ai(0), as a function of an observable, O, in a given bin, i, is 
Ai(0) = Nl ec {0)/Nf n {0), where Nf n (0) is the number of generated MC events and 
Nl ec iO) is the number of reconstructed events passing all the selection requirements. 

Differential cross sections as a function of O in a given bin i were obtained using the 
expression 

d(7i Ni 

dO = B C Ai{oy 

where N is the number of signal events, reconstructed in each bin after subtraction of the 
estimated contribution from the diffractive proton-dissociative events, B the branching 
ratio (5.88 ± 0.10)% [35] and C the integrated luminosity. 

The background from if>' photoproduction is expected to be 15% [49]; this expectation 
was confirmed by a direct measurement of the ip' to J /if) cross section ratio [13]. Re- 
stricting the phase-space region in this analysis similar to that for photoproduction, 
50 < W < 180 GeV and 0.55 < z < 0.9, the number of observed if)' events was con- 
sistent with the expectation from the if)' to J/if> ratio measured in the photoproduction 
regime. The contribution of J / if) mesons from if)' decays was assumed to yield the same 
kinematic distributions as the dominant direct J /if) contribution and, therefore, the the- 
oretical predictions for J /if) production were scaled up by 15%. This change is small 
compared to the normalisation error of the LO NRQCD predictions. 

Monte Carlo studies showed that the contribution from 5-meson decays into J/ if) was 
concentrated at low— z values and small elsewhere. For 0.1 < z < 0.4, this contribution 
can be as large as 20%. The beauty contribution was estimated using the Rapgap MC 
and added to the J/if> predictions. This change is small compared to the normalisation 
uncertainty of the LO NRQCD predictions. 

The J /if) meson can be produced via Xc radiative decays, Xc — > J/4>1- While Xc mesons 
can be produced copiously in hadron-hadron collisions through gg, gq and qq interactions, 
Xc production via photon-gluon fusion is forbidden at LO in the CS model. This leaves 
only resolved photon processes, strongly suppressed at non-zero photon virtuality, or CO 
processes as sources of Xc production. However, the ratio of the Xc to J / if) from the CO 
processes is expected to be below 1% [50]. This contribution was therefore neglected. 

The effect of the LO electroweak corrections was studied using the Heracles [51] MC 
program. The open charm DIS cross section was evaluated using the Rapgap [45] MC 
program with and without radiative corrections, as calculated by Heracles, in a W-Q 2 
grid. The measured cross sections were then corrected to the QED Born level using the 
Heracles predictions. In the region covered by the data, this correction was —2% on 
average and always below 7% in absolute value. 
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6 Systematic uncertainties 



The systematic uncertainties of the measured differential cross sections were determined 
by changing the selection cuts or the analysis procedure in turn and repeating the ex- 
traction of the differential cross sections. The resulting uncertainty on the total cross 
section is given in parentheses. The following categories of systematic uncertainties were 
considered: 

• scattered electron reconstruction: these uncertainties were evaluated as described else- 
where [52] (2%); 

• CAL energy scale and resolution simulation: these uncertainties were evaluated as 
described elsewhere [52] (2%); 

• tracking: the resolutions on track momenta and angles were varied by ± 20% of their 
values and the magnetic field by ± 0.3% (1%); 

• muon reconstruction: the uncertainty of the muon acceptance, including those of 
the efficiency of the muon chambers, the trigger selection algorithms and the offline 
reconstruction, was obtained from a study based on an independent dimuon sample 
at high Q 2 , performed following the method discussed elsewhere [53] (6%); 

• fitting procedure: the invariant-mass range and the functional form of the background 
were varied (2%); 

• simulation of the process 7*g — > J/ijjg: the Cascade MC rather than the Epjpsi MC 
was used to calculate acceptances (5%); 

• subtraction of the remaining diffractive proton-dissociative admixture: Diffvm rather 
than Epsoft was used to perform the subtraction of the proton diffractive events (3%). 

These estimations were also made in each bin of the differential cross sections. All of the 
above individual sources of systematic uncertainty were added in quadrature. 

The following sources resulted in an overall shift of the cross section and were therefore 
treated as normalisation uncertainties: 

• the integrated luminosity determination has an uncertainty of 2%; 

• the branching ratio of J/ip — > has an uncertainty of 1.7% [35]. 

The normalisation uncertainties were not included in the total systematic uncertainty. 

7 Results 

The cross section for the process ep —* eJ/ipX in the kinematic region 2 < Q 2 < 80 GeV 2 , 
50 < W < 250 GeV, 0.2 < z < 0.9 and -1.6 < F lab < 1.3 is 
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302 ± 23 (stat.) 



+28 
-20 



(syst.) pb, 



where the first uncertainty is statistical and the second systematic. In Figs. El ID and 
the differential cross sections as a function of z, Q 2 , W, p^, Y*, logMj^ and the rapidity 
of the hadronic system X, Yx, are shown. They are compared to the predictions of a 
NRQCD model [54], a CS model with k T factorisation (LZ) [18] and to the Cascade 
MC. The beauty contribution, estimated using the Rapgap MC, is also shown separately 
in Figs. Efc and 0^. All differential cross sections and normalised cross sections are given 
in Tables □ and H 

The uncertainties for the CS and CO NRQCD predictions correspond to variations of 
the charm-quark mass (m c = 1.5 ±0.1 GeV) and of the renormalisation and factorisation 
scales from V2V Q 2 + M 2 to 2V Q 2 + M|. The uncertainty on the long-distance matrix 
elements and the effect of different choices of parton distribution functions (default set 
is MRST98LO) are also taken into account. The bands in the figures shows all these 
uncertainties added in quadrature. 

In general, the CSM is consistent with the data. The predictions including both CS 
and CO contributions are higher than the data, especially at high z and low p^. At 
high values of the agreement with the data is reasonable. The prediction does not 
describe the shapes of the z, Y* , logMj^ and Yx distributions. Previous photoproduction 
results [12, 13] showed that the agreement between data and theory at high z can be 
improved using resummed LO NRQCD predictions [55]. It should be noted that, in 
photoproduction, inclusion of the NLO corrections to the CSM, not available for DIS, 
significantly improved the description of the data. 

For the LZ ^-factorisation predictions, the parametrisation, KMS [56], of the uninte- 
grated gluon density was used. The charm-quark mass was set to m c = 1.4 GeV, which is 
the mass used in the KMS parametrisation. The renormalisation and factorisation scales 
were both set to p = kr for kr > 1 GeV. For kr < 1 GeV the scales were fixed at 1 GeV. 
Calculations based on the /c^-factorisation approach give a reasonable description of the 
data both in shape and normalisation. 

The data are also compared with the predictions of the Cascade MC using the kr- 
factorisation approach, where gluons are treated according to the CCFM evolution equa- 
tions. These predictions were obtained by setting the charm-quark mass to 1.5 GeV, the 



evolution scale of the strong coupling constant to the J/ip transverse mass, y M 2 
and using the unintegrated gluon-density parametrisation " J2003 set 2" . The CASCADE 
MC is above the data for z > 0.45 and for W < 175 GeV. 

In order to compare the present measurements directly to the HI results [14], differential 
cross sections were determined in the kinematic range 2 < Q 2 < 100 GeV 2 , 50 < W < 
225 GeV, 0.3 < z < 0.9 and p^ > 1 GeV 2 ; all ZEUS differential cross sections and 
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normalised cross sections are given in Table El The results of this comparison are shown 
in Fig. The present results are in agreement with those from HI. In Fig. the 
ZEUS data are in better agreement with the CSM prediction than in Fig. This is a 
consequence of the p^ > 1 GeV 2 cut used in Fig. combined with the fact that the CS 
prediction underestimate the data at high p^, as seen in Fig. 0^. 



8 Conclusions 



Inelastic J/ip production in DIS has been measured in the kinematic region 2 < Q 2 < 
80 GeV 2 , 50 < W < 250 GeV, 0.2 < z < 0.9 and -1.6 < Y lab < 1.3. The data are in agree- 
ment with the HI results in the kinematic region 2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, 
0.3 < z < 0.9 and p^ 2 > 1 GeV 2 . The data are compared with LO NRQCD predictions, 
including both CS and CO contributions, and fey-factorisation calculations. Calculations 
of the CS process generally agree with the data, whereas inclusion of CO terms spoils this 
agreement. 
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z range 


da/dz (pb) 


1/ada/dz 


0.20-0.45 
0.45-0.60 
0.60-0.75 
0.75-0.90 


309 ± 611fJ 
428 ± 62+^ 
568 ± 651^ 
526 ± 66l$ 


i.oi ±0.1618:31 
i.4o ±o.i9i8:™ 

1.86 ± 0.20l8:?3 
1.72 ± 0.20l8:ls 


W range (GeV) 


da/dW (pb/GeV) 


1/ada/dW 


50-100 
100-125 
125-175 
175-250 


1.73 ± 0.25±g;fg 
2.44 ± 0.32+°; 2 ^ 
1.43±0.20lg;J| 
1.17 ±0.22+°;^ 


0.0056 ± 0.000718:8805 

0.0080 ±0.0010l8:888s 

0.0047 ±0.000618:8883 
o.oo38 ± o.ooo6i8:8881 


Q 2 range (GeV 2 ) 


da/dQ 2 (pb/GeV 2 ) 


1/ada/dQ 2 


2-4 
4-8 
8-16 
16-80 


66.9 ± 8At 7 6 i 
18.3 ± 2.7±i;| 
6.3 ± l.Otol 
0.66 ±0.12+°;^ 


0.223 ± 0.019l8:8?2 

0.0609 ± o.oo79i8:88l 

_ „ _ . , „ „ „^„_L(1 (1(111 

0.0211 ± 0.00 3218:88^ 
0.00221 ± o.ooo38i8:888i4 


range (GeV 2 ) 


da/dpf (pb/GeV 2 ) 


l/(jdcr/dp^ 2 


0- 1 

1- 5 
5-16 

16-100 


80 ± 14ljj 
40.1 ±4.ll^ 
3.81 ± 0.70l{j;^ 
0.280 ± 0.0511^ 


0.269 ± 0.04ll8:83 2 
0.1345 ±0.0096l8:88i4 
0.0128 ±0.0022l8:88o8 
0.00094 ± 0.00017l8:8888g 


Y* range 


da/dY* (pb) 


1/ada/dY* 


1.75-2.60 
2.60-3.00 
3.00-3.40 
3.40-4.00 


80 ± 16±? 
212 ± 28l?8 
211 ± 25+l| 

94 ± Ut l 9 8 


0.274 ± 0.04518:821 

0.722 ± 0.08318:831 

0.716 ± 0.07718:81 
0.321 ± 0.04518:821 



Table 1: Differential cross sections and normalised differential cross sections in 
the kinematic region 2 < Q 2 < 80 GeV 2 , 50 < W < 250 GeV, 0.2 < z < 0.9 and 
— 1.6 < Yi a b < 1.3 as a function of z, W, Q 2 , p^ 2 and Y* . The first uncertainty is 
statistical and the second is systematic. Overall normalisation uncertainties due to 
the luminosity measurement (±2%) and to the J/ip decay branching ratio (1.7%) 
are not included in the systematic error. 
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log(Mj./GeV 2 ) range 


da/d\og(M x /GeV 2 ) (pb) 


l/ada/d\og(M x /GeV 2 ) 


3.00-3.55 
3.55-3.85 
3.85-4.10 
4.10-4.50 


156 ± 18+^ 
208 ± 27± 2 ^j 
270 ± 40±|f 
164 ± 31±?| 


0.556 ± 0.057 + °;° 1 ^ 
0.740 ± 0.0911^ 

0.96 ±0.131^ 
0.581 ± 0.092+°;°^ 


Yx range 


da/dYx (pb) 


1/ada/dYx 


2.20-2.78 
2.78-3.05 
3.05-3.37 
3.37-4.05 


112 ±21±\l 
243 ± 37j|jj 
203 ± 26±f 5 
143 ± 16t^ 


0.383 ± 0.061+°;°| 

0.83 ±0.1lj£$ 
0.692 ± 0.083t{j;g^ 
0.488 ± 0.047±g;gH 



Table 2: Differential cross sections and normalised differential cross sections in 
the kinematic region 2 < Q 2 < 80 GeV 2 , 50 < W < 250 GeV , 0.2 < z < 0.9 and 
— 1.6 < Yi ab < 1.3 as a function of\og(M x / GeV 2 ) and Y x . The first uncertainty is 
statistical and the second is systematic. Overall normalisation uncertainties due to 
the luminosity measurement (±2%) and to the J/ip decay branching ratio (1.7%) 
are not included in the systematic error. 
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z range 


da/dz (pb) 

/ \ J- / 


1/adcr/dz 


0.30-0.45 
0.45-0.60 
0.60-0.75 
0.75-0.90 


246 ± 60+| 

317 ±50j|2 
430 ± 56l|J 
392 ± 57lg 


1.18 ±0.25!°;^ 
1.53±0.22l°;{3 

2.07±0.23lHi 
1.89 ± 0.24l°; 22 


Pt range (GeV 2 ) 


do /dp* 2 (pb/GeV 2 ) 


l/odcr/dp*? 


1-5 
5-16 
16-40 


36.4 ±3.71^ 
3.65 ± 0.711°;^ 
0.92±0.18l^ 


0.1752 ±0.0092lg;gg^ 
0.0176 ± 0.0030l^ 22 
0.00443 ± 0.00083lg;E 


Y* range 


da/dY* (pb) 


1/ada/dY* 


2.00-2.60 
2.60-3.00 
3.00-3.40 
3.40-4.00 


66 ± 14+ r 10 
137 ± 20l^ 
144 ± 19l{? 
61 ± 14l^ 2 


0.351 ± 0.066l°;° 2 ^ 
0.728 ± 0.0941^2 
0.762 ± 0.09118;^ 
0.323 ± 0.0641°;°^ 



Table 3: Differential cross sections and normalised differential cross sections in 
the kinematic region 2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV, 0.3 < z < 0.9 and 
Pt > 1 GeV 2 as a function of z, p*? and Y* . The first uncertainty is statistical and 
the second is systematic. Overall normalisation uncertainties due to the luminosity 
measurement (±2%) and to the J/ip decay branching ratio (1.7%) are not included 
in the systematic error. 
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Q 2 < 80 GeV 2 , 50 < W < 250 GeV, 0.2 < z < 0.9 and -1.6 < Y lab < 1.3 plotted 
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Figure 3: Differential cross sections for the reaction ep — > eJ/ipX in the 
kinematic region 2 < Q 2 < 80 GeV 2 , 50 < W < 250 GeV, 0.2 < z < 0.9 and 
— 1.6 < Yi a b < 1.3 as a function of (a) z, (c) Q 2 and (e) W . The inner error bars 
of the data points show the statistical uncertainty; the outer bars show statistical 
and systematic uncertainties added in quadrature. The data are compared to LO 
NRQCD predictions, a LO CS calculation, a prediction in the kx -factorisation 
approach within the CSM and the Cascade MC predictions, (b), (d) and (f) show 
the data and the theoretical predictions normalised to unit area. 
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Figure 4: Differential cross sections for the reaction ep — > eJ/ipX in the 
kinematic region 2 < Q 2 < 80 GeV 2 , 50 < W < 250 GeV , 0.2 < z < 0.9 and 
— 1.6 < Yiah < 1.3 as a function of (a) and (c) Y* . The inner error bars 
of the data points show the statistical uncertainty; the outer bars show statistical 
and systematic uncertainties added in quadrature. The data are compared to LO 
NRQCD predictions, a LO CS calculation, a prediction in the kj '-factorisation 
approach within the CSM and the Cascade MC predictions, (b) and (d) show the 
data and the theoretical predictions normalised to unit area. 
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Figure 5: Differential cross sections for the reaction ep — > eJ/ipX in the 
kinematic region 2 < Q 2 < 80 GeV 2 , 50 < W < 250 GeV , 0.2 < z < 0.9 and 
— 1.6 < Ffa6 < 1.3 as a function of (a) log Mj^ and (c) Y x . The inner error bars 
of the data points show the statistical uncertainty; the outer bars show statistical 
and systematic uncertainties added in quadrature. The data are compared to LO 
NRQCD predictions, a LO CS calculation, a prediction in the kx -factorisation 
approach within the CSM and the Cascade MC predictions, (b) and (d) show the 
data and the theoretical predictions normalised to unit area. 
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Figure 6: Differential cross sections for the reaction ep — > eJ/ipX in the 
kinematic region 2 < Q 2 < 100 GeV 2 , 50 < W < 225 GeV , 0.3 < z < 0.9 and 
p^ > 1 GeV 2 as a function of (a) z, (c) p*^ and (e) Y* . The inner error bars 
of the data points show the statistical uncertainty; the outer bars show statistical 
and systematic uncertainties added in quadrature. The ZEUS and HI data are 
compared to LO NRQCD predictions, a LO CS calculation, a prediction in the fcy- 
factorisation approach within the GSM and the Cascade MC predictions. The HI 
data points are plotted at the mean value of the data in each interval [14]- The 
ZEUS data for the p^ differential cross section are plotted at the weighted mean, 
for each bin, of the Epjpsi MC prediction, (b), (d) and (f) show the data and the 
theoretical predictions normalised to unit area. 
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